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Abstract 

The effects of delayed feedback terms on nonlinear oscillators has been extensively 
studied, and have important applications in many areas of science and engineering. 
We study a particular class of second-order delay-differential equations near a point of 
triple-zero nilpotent bifurcation. Using center manifold and normal form reduction, we 
show that the three-dimensional nonlinear normal form for the triple-zero bifurcation 
can be fully realized at any given order for appropriate choices of nonlinearities in the 
original delay-differential equation. 
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1 Introduction 



Delay-differential equations are used as models in many areas of science, engineering, eco- 
nomics and beyond [21 HH HH1 [13 UHl HH1 [201 EH [23l [21]. It is now well understood that 
retarded functional differential equations (RFDEs), a class which contains delay-differential 
equations, behave for the most part like ordinary differential equations on appropriate 
infinite-dimensional function spaces. As such, many of the techniques and theoretical re- 
sults of finite-dimensional dynamical systems have counterparts in the theory of RFDEs. In 
particular, versions of the stable/unstable and center manifold theorems in neighborhoods of 
an equilibrium point exist for RFDEs [13J. Also, techniques for simplifying vector fields via 
center manifold and normal form reductions have been adapted to the study of bifurcations 
in RFDEs 0(10]. 

One of the challenges of applying these finite-dimensional techniques to RFDEs lies in the 
so-called realizability problem. This problem stems from the fact that the procedure to reduce 
an RFDE to a center manifold often leads to algebraic restrictions on the nonlinear terms in 
the center manifold equations. Specifically, suppose B is an arbitrary mxm matrix. For the 
sake of simplicity, suppose additionally that all eigenvalues of B are simple. Let C([—r, 0], R) 
be the space of continuous functions from the interval [— r, 0] into R, and for any continuous 
function z, define z t e C([-r, 0],R) as z t (6) = zii + 6), -r < 9 < 0. It is then possible [11] 
to construct a bounded linear operator L : C([— r, 0],R) — > R such that the infinitesimal 
generator A for the flow associated with the functional differential equation 

z(t)=Lzt (1.1) 

has a spectrum which contains the eigenvalues of B as a subset. Thus, there exists an Tri- 
dimensional subspace P of C([— r, 0],R) which is invariant for the flow generated by A, and 
the flow on P is given by the linear ordinary differential equation (ODE) 

x = Bx. 

Now, suppose fll.ip is modified by the addition of a nonlinear delayed term 

z(t) = L z t + az{t - t) 2 , (1.2) 

where a 6 IS is some coefficient and r G [0, r] is the delay time. Then the center manifold 
theorem for RFDEs [32] can be used to show that the flow for (11.21) admits an m-dimensional 
locally invariant center manifold on which the dynamics associated with (11. 2p are given by 
a vector field which, to quadratic order, is of the form 

x = Bx + ag(x), (1.3) 
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where g : M. m — > M. m is a fixed homogeneous quadratic polynomial which is completely 
determined by L and r, and a is the same coefficient which appears in (II. 2p . We immediately 
notice that for fixed £ and r, (jl.3p has at most one degree of freedom in the quadratic term, 
corresponding to the one degree of freedom in the quadratic term in (|1.2j) . However, whereas 
one degree of freedom is sufficient to describe the general scalar quadratic term involving 
one delay in (11.21) . it is largely insufficient (if m > 1) to describe the general homogeneous 
quadratic polynomial / : K' m — > M. m . Therefore, there exist m-dimensional vector fields 
x = Bx + f(x) (where / is homogeneous quadratic) which can not be realized by center 
manifold reduction (II. 3p of any RFDE of the form (II. 2p . The realizability problem has 
received considerable attention in the literature [3j El [6j [TTJ [12] . 

In this paper, we will be interested in a realizability problem for a class of second-order 
scalar delay-differential equations of the form 

x(t) + bx(t) + ax(t) - F(x(t),x(t)) = ax(t - r) + I3x{t - r) + G(x(t - r),x(t - r)), (1.4) 

where a, b, a and f3 are real parameters, r > is a delay term, and the nonlinear functions 
F and G are smooth and vanish at the origin, along with their first order partial derivatives. 
This class contains many interesting applications which have been studied in the literature, 
including Van der Pol's oscillator with delayed feedback [TJ [7J [151 1221 [23], as well as models 
for stabilization of an inverted pendulum via delayed feedback [IS] . 

Both the Van der Pol oscillator [23] and the inverted pendulum system [IS] have been 
shown to possess points in parameter space where a bifurcation via a non-semisimple triple- 
zero eigenvalue occurs. In [IS], this bifurcation is in fact characterized as the organizing cen- 
ter for their model, since it includes in its unfolding Bogdanov-Takens and steady-state/Hopf 
mode interactions and the associated complex dynamics of these codimension two singulari- 
ties. As far as we are aware, a complete theoretical analysis and classification of all possible 
dynamics near the non-semisimple triple-zero bifurcation has yet to be done, although a 
rather thorough investigation was undertaken in [8j. Numerical tools are used in [IS] to il- 
lustrate the complexity of this singularity in their model, including many global bifurcations. 
It is stated in [IS] that because of the presence of invariant tori, a full versal unfolding of 
the triple-zero singularity must include terms other than those appearing at cubic order in 
their model, and conclude by wondering whether full realizability of the nonlinear normal 
form for the triple-zero bifurcation is possible for their delay-differential equation. 

Other relevant work includes [I], where the authors study a class of coupled first-order 
delay-differential equations which includes (II. 4p as a special case (if one writes (II. 4p as a first 
order system), and compute quadratic and cubic normal form coefficients in term of DDE 
coefficients for both non-semisimple double-zero and triple-zero bifurcations. Higer-order 
terms for these normal form are not considered. 
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In this paper, we will first show that the non-semisimple triple-zero singularity occurs 
generically in ( II. 4p , and then prove that the full nonlinear normal form for the non-semisimple 
triple-zero bifurcation, at any prescribed order, can be realized by center manifold normal 
form reduction of ( II ,4p for appropriate choices of nonlinear functions F and G. In section 2, 
we present the functional analytic framework in which we will study this problem. Section 3 
gives a brief summary of the center manifold and normal form procedure for RFDEs which 
was developed by Faria and Magalhaes [HI HD]. Our main result is stated and proved in 
section 4. We end with some concluding remarks in section 5. 

2 Functional analytic setup 

As mentioned in the introduction, we consider a general class of second order nonlinear 
differential equations for the real- valued function x(t) of the form (II. 4p . which we rewrite as 
a first order system 



x(t) = y(t) 

y(t) = -ax(t) + ax(t-r) - by(t) + /3 y(t - r) + F(x(t),y(t)) + G(x(t - r),y(t - r)). 



where a, 6, a and (3 are real parameters, r > is a delay term, and the nonlinear functions 
F and G are smooth and vanish at the origin, along with their first order partial derivatives. 
In many applications, we have a > 0, so we will assume this condition throughout (although 
other cases of a can be treated in a similar manner). 

The characteristic equation corresponding to (12. ip is P(X) = 0, where 



(2.1) 



P(A) = A 2 + bX + a - (a + (3X)e 



Ar 



(2.2) 



A straightforward computation shows that when 



a = a 



r = to = - + 
a 



{3 ^(3 2 + 2a 



a 



(2.3) 



b = f3 — ar 



3(3 2 7^ 2a (or equivalently ar ^ 3/3) 
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then -P(O) = P'(O) = P"(0) = 0, P"'(0) 7^ 0, and P has no other roots on the imaginary axis. 
Therefore is a triple eigenvalue for the linearization of (12. ip at the origin, with geometric 
multiplicity one. 

For the parameter values (12.31) . we write (12. ip as 

x(t) = y{t) 

y(t) = -a(x(t)-x(t-r )) + ar o y(t)-/3(y(t)-y(t-r )) (2.4) 

+F(x(t),y(t)) + G{x{t - r ),y{t - r )). 

Let C = C ([— ro,0] , M 2 ) be the Banach space of continuous functions from [— To,0] into 
R 2 with supremum norm. We define z t G C as 

z t (6) = z(t + 9)= + ),-t o <0<O. 

We view (I2.4p as a retarded functional differential equation of the form 

z{t) = JCz t + 3 : (zt), (2.5) 
where L : C M 2 is the bounded linear operator 

' 1 \ , , / 



and IF is the smooth nonlinear function from C into IR 2 

^ = ( P(0(O)) + G(0(-r o )) 

Let A be the infinitesimal generator of the flow for the linear system z = L Zt, with 
spectrum a(A) D {0}, and P be the three-dimensional invariant subspace for A associated 
with the eigenvalue 0. Then it follows that the columns of the matrix 

*=( 1 9 ¥ 2 

\0 1 9 
form a basis for P. 

In a similar manner, we can define an invariant space, P*, to be the generalized eigenspace 
of the transposed system, A T associated with the triple nilpotency having as basis the rows 
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of the matrix \1/ =col(ipi, . . . ,ip m ). Note that the transposed system, A T is defined over a 
dual space C* = C ([0, r ] , R 2 ) , and each element of \1/ is included in C*. The bilinear form 
between C* and C is defined as 

o 8 

(lM) = V (0)0(0)- / U{C-e) [dr,(9)\ 0(C) dC- (2-6) 



-r 

Note that $ and \1> satisfy $ = 5$, \i> = —^B, where S is the 3 x 3 matrix 



5=001. (2.7) 




We can normalize \I/ such that (\l/, $) = /, and we can decompose the space C using the 
splitting C = P © Q, where the complementary space Q is also invariant for A. 

Faria and Magalhaes 0, [TO] show that f)2.5p can be written as an infinite dimensional 
ordinary differential equation on the Banach space BC of functions from [— tq,0] into M 2 
which are uniformly continuous on [— r ,0) and with a jump discontinuity at 0, using a 
procedure that we will now outline. Define Xq to be the function 



Xo(0) 



1 

1 







9 = 

-t <6 < 0, 



then the elements of BC can be written as £ = tp + X X, with 99 6 C and A G M 2 , so that 
BC is identified with C x R 2 . 

Let 7r : .BC — )■ P denote the projection 

7r(p + X A) + tt(0)A], 

where tp E C and A G R 2 . We now decompose z$ in (I2.5P according to the splitting 

BC = P® ker tt, 

with the property that Q C ker7r, and get the following infinite-dimensional ODE system 
which is equivalent to ( 12. 5ft : 

u = Bu + ^(0) + 

d (2-8) 
— v = A Q iv+(I -7r)X 3($ u + v), 

(3 



where u G M 3 , v E Q 1 = Q H C 1 , (C 1 is the subset of C consisting of continuously differen- 
tiable functions), and Aqi is the operator from Q l into ker7r defined by 

A Q np = ip + X [L f - (f(0)]. 

3 Faria and Magalhaes normal form 

Consider the formal Taylor expansion of the nonlinear terms 7 in (12. 5p 

j(0) = ]T^(0), 0g c, 

i>2 

where 3^(0) = Vj(<fi, . . . ,<fi), with Vj belonging to the space of continuous multilinear sym- 
metric maps from C x - • • x C ( j times) to R 2 . If we denote fj = (fj, ff), where 



fj(u,v) = VMfjiQu + v) 
ff(u,v) = (I-n)X ^u + v), 



then (12. 8p can be written as 



a = Bu + yj fj{u, r 



A Q iv + }_^ ff(u,v). 
i>2 



~r v 
dt 



It is easy to see that the non-resonance condition of Faria and Magalhaes (Definition 
(2.15) of [TO] ) holds. Consequently, using successively at each order j a near identity change 
of variables of the form 

(u,v) = (u,v) + U j (u) = (u,v) + {Uj(u),Uf(u)), (3.2) 

(where £7' 2 are homogeneous degree j polynomials in the indicated variable, with coefficients 
respectively in M 3 and Q l ) system (13. ip can be put into formal normal form 

u = Bu + ^2gj(u,v) 

d A ^ 2, , ( 3 ' 3 ) 

— v = A Q iv + ^ gj (u,v) 
i>2 
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such that the center manifold is locally given by v = and the local flow of (12. 5 j) on this 
center manifold is given by 

u = Bu + ^2g^(u,0). (3.4) 

i>2 

The nonlinear terms in (I3.4p are in normal form in the classical sense with respect to the 
matrix B. 



4 Realizability of the normal form for the triple-zero 
nilpotency 

It was shown in [SJ [21] that the classical normal form for the general nonlinear vector field 

«i \ /010\/ui\ / r 1 (u 1 ,u 2 ,u 3 ) 
u 2 I = 1 \ I u 2 + r 2 (u 1 ,u 2 ,u 3 ) 
u 3 J \ / \ ii 3 / \ r 3 («i,tt2,w 3 ) 

dv ■ 

(where r,-(0, 0, 0) = 0, 7^(0, 0, 0) = 0, for j, k = 1, 2, 3) is 

Wi = u 2 
u 2 = u 3 

j>2 \ i=0 i=0 / 

where 

N = 1, J = |(j — 1), / = «/, when j is odd, 

N = 2, J = 2 — 1, / = J + 1, when j is even. 

Thus, if 5 is the matrix (12. 7p . and if- 7 is the space of homogeneous polynomial mappings 
of degree j > 2 from K 3 into M 3 , then the homological operator Lb = Dh(u) ■ Bu — B ■ h(u) 
acting on W is such that 

H j = L B {H j )®Wj, (4.1) 
where Wj C iP is the subspace of dimension 

3j + 2 



2 

3j + 3 



if j is even 
if j is odd 



spanned by 




j—i i 

U\ u 2 



i = 0,...,3 } |J 






j—l—i i+1 



z = 0, . . . , |(j — 2) j even 

» = o, . . . , |0" - 1) ^' odd - 



Now, if F(zi,z 2 ) and Cr(zi,Z2) are real-valued functions such that F, G and their first- 
order partial derivatives vanish at the origin, then we may write the Taylor series 

F(z h z 2 ) = ^2 F i{ z ^ z ?)i G(z 1 ,z 2 ) = ^2 Gj(z u z 2 ), 

j>2 j>2 

where the Fj and Gj are homogeneous degree j polynomials. The first equation in (13. ip then 
reduces to 



u 



i>2 



/ o 

«i (Fj((ui, u 2 ) + v(0)) + Gj((ui - t u 2 + |t 2 m 3 , u 2 - t u 3 ) + u(-r ))) 
\ k 2 {Fj((ui, u 2 ) + v(0)) + Gj{{ui - r u 2 + |t 2 m 3 , u 2 - r u 3 ) + u(-r ))) / 



(4.2) 



where 



3(ar - 4/3) 



2r (ar -3/3) 2 ' 
We note that when v = in (14.21) . then 

Fj{ui,u 2 ) + Cr,-(ui - r u 2 + |r 2 M 3 , w 2 - t u 3 ) 



«2 



r 2 (ar - 3/9) 



^0. 



A/(l*l, U 2 ) + «lW3%-2(Ml, U 3 ) + 
u 2 u 3 ej- 2 (ui, U 2 , U 3 ) + W 2 D j _ 2 (M 3 ), 



where 



Aj (u u u 2 ) = F j (u 1 ,u 2 ) + Gjiv,! - t u 2 ,u 2 ) 



(4.3) 



u 1 u 3 'B j _ 2 (u 1 , us) = &j(ui + ^r 2 w 3 , -t u 3 ) - Gj(ut, 0) - G^-TqU^ -r u 3 ) (4.4) 

«2%e j _ 2 (Ml, U 2 ,U 3 ) = Gj(Ui - T U 2 + ^TqU 3 ,U 2 - TqU 3 ) - Gj(ui - T U 2 ,U 2 )- 

GjiUx + ~TqU 3 , -T Q U 3 ) + Gj(lli, 0) 

ulDj- 2 (u 3 ) = Gj(]-T$U 3 , -t q u 3 ). 
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Lemma 4.1 For a given integer j > 2, let ((ui,u 3 ) be a homogeneous degree j polynomial 
such that £(0,113) — C( w i)0) = 0- Then there exists a homogeneous polynomial of degree j, 
£(iti, u 3 ) such that 

I 1 

((ux, u 3 ) = £(ui + -TqU 3 , -t u 3 ) - £(«i, 0) - £(-t 2 « 3 , -r u 3 ). (4.5) 

Proof of lemma: If we write 

3 

£(ui,u 3 ) = lj-i,i u l 



■1 % 
u 3 



i=0 



then a lengthy but straightforward computation shows that 
+ ^r 2 u 3 , -t u 3 ) - 0) - £(^r 2 u 3 , -r u 3 ) = 

'i-l / ■ \ /1 \ k \ i-i /j-i-l 



->» (E ( j ) (r») «!-*<*) + E f E ( J * ' ) (r») Htf**V 



where 



fc J k\(j-k)l 
Now, since £(0,m 3 ) = ((ui,0) = 0, we have that 

j'-i 

C(«l,%) = ^ ej-i,iM{~ Ug. 

1=1 

Thus, we see for example that we may solve (14.51) by arbitrarily setting 7^0 = 0, 7oj = 0, 
and choosing 7j-i,i, i = 1, . . . , j — 1 such that the following triangular linear algebraic system 
is satisfied 

3 1 ) (~ r o) 7j-i,i = 

( 3 1 1 ) (r° 2 ) ( ~ ro) 7j '- 1 - 1 + ( J 2 ) (_r ° )2 7j - 2 ' 2 = e '- 2 ' 2 



-Tn 



(-To) 7i-l,l + 



+ 
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This ends the proof of the lemma. I 

Now, recalling the splitting (14. ip . let Qj(u\, u 2 , u 3 ) be a homogeneous degree j polynomial 
such that Qj £ Wj. We may write 









^l%Sj-20l,%) 



where qj is a homogeneous polynomial of degree j, and Sj- 2 is a homogeneous degree j — 2 
polynomial. 

The degree j term in (14.21) for v = can be written as 



/ 







\ 



Kl (F j (ui,U 2 ) + Gj(Ui - T U 2 + |To« 3 , W 2 - T U 3 )) 



\ K 2 «2) + Gj(ui - T U 2 + ^TqU 3 , U 2 - T U 3 )) ) 

( \ 

Rj(ui,u 2 ,u 3 ) + I 



+ 








K 2 Aj(u\, u 2 ) 



K 2 UxU 3 'B j -2(ui,U 3 ) 



where Rj is in the range of the homological operator, Rj C Le(i/ J ), and Aj and "Bj- 2 are 
as in ( I4.3p and ( I4.4p . Using the previous lemma, we know that if we choose Gj such that 
k 2 "Bj-2 = Sj_ 2 , and then set k 2 Fj(ui, u 2 ) = qj(ux,u 2 ) — k 2 Gj(ui — t u 2 , u 2 ), then 













(4.6) 



K, 2 Aj(ui,U 2 ) ) \ K 2 UlU 3 'Bj-2(Ul,U 3 ) 

We can now state and prove the following realizability theorem: 
Theorem 4.2 Given an integer £ > 2 and a polynomial vector field on M. 3 of the form 



u 



Bu + Wj(u), 



(4.7) 



i=2 



where B is the matrix (2.1) and Wj £ Wj as in j4-l\) , there exist polynomial functions F and 
G in ( fjO ) such that the Faria and Magalhaes center manifold and normal form reduction 
(ggP 0/ (fg up to order £ is $7fy . 
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Proof: Applying successively at each order j (from j = 2 to j = £) a near identity change 
of variables of the form (I3.2p . and setting v = 0, we transform (14 .2p into 



u 



3=2 






kqAAui, u 2 ) 






K 2 U 1 U 3 'B i ^ 2 (u 1 ,U 3 ) 



+ A j (u 1 ,u 2 ,u 3 ) 



+ 0{\u 



(4.8) 

where A 2 (ui,u 2 ,u 3 ) = and for j > 3, Aj G Wj is an extra contribution to the terms of 
order j coming from the transformation of the lower order (< j) terms. Therefore, we set 

Qj(u) = Wj(u) - Aj(u), j = 2,...,£, 

and use (14. 6 j) to conclude that the truncation of (14. 8 p at order £ is (14. 7p . ■ 



5 Conclusion 

In this paper, we have solved the realizability problem for the normal form of the non- 
semisimple triple- zero singularity in a class of delay differential equations ( jl .41) which includes 
delayed Van der Pol oscillators, as well as certain models for the control of an inverted 
pendulum as special cases. It is apparent from the complexity of the dynamics of ( 11.411 
near the triple-zero nilpotency reported in previous work [U El [18] that high-order normal 
forms will be required for a complete classification of this singularity. Although such a 
complete classification of the dynamics near a triple-zero nilpotency is beyond the scope of 
this paper, our results allow us to conclude that the full range of complexity of this singularity 
is attainable within the class of delay-differential equations (11. 4p . 

Although we have not done so, we believe that the results in this paper could be suitably 
generalized to studying realizability of higher order nilpotencies in higher-order scalar delay- 
differential equations such as 

ra—l n— 1 

x (n \t)+J2 a j x®(t)-F(x(t),...,x<> n - 1 ')(t)) = ^ OL j x^{t-r)+G{x{t-T),...,x { - n - 1 \t-T)) 

j=0 j=0 

where n > 3 is an integer. 

Acknowledgments 

This research is partly supported by the Natural Sciences and Engineering Research 
Council of Canada in the form of a Discovery Grant. 



12 



References 



[1] F.M. Atay. Van der Pol's oscillator under delayed feedback. J. Sound and Vibration 
218, (1998) 333-339. 

[2] A. Beuter, J. Belair and C. Labrie. Feedback and delays in neurological diseases : a 
modeling study using dynamical systems. Bulletin Math. Biology 55, (1993) 525-541. 

[3] P.-L. Buono and J. Belair. Restrictions and unfolding of double Hopf bifurcation in 
functional differential equations. J. Differential Equations 189, (2003) 234-266. 

[4] S.A. Campbell and Y. Yuan. Zero singularities of codimension two and three in delay 
differential equations. Nonlinearity 21, (2008) 2671-2691. 

[5] Y. Choi and V.G. LeBlanc. Toroidal normal forms for bifurcations in retarded functional 
differential equations I: multiple Hopf and transcritical/multiple Hopf interaction. J. 
Differential Equations 227, (2006) 166-203. 

[6] Y. Choi and V.G. LeBlanc. Toroidal normal forms for bifurcations in retarded func- 
tional differential equations II: saddle- node/multiple Hopf interaction. Dynamics of 
Continuous, Discrete and Impulsive Systems A 15 (2008), 251-276. 

[7] J.C.F. de Oliveira. Oscillations in a van der Pol equation with delayed argument. J. 
Math. Anal. Appl. 275, (2002) 789-803. 

[8] F. Dumortier and S. Ibahez. Nilpotent singularities in generic 4-parameter families of 
3-dimensional vector fields. J. Differential Equations 127, (1996) 590-647. 

[9] T. Faria and L.T. Magalhaes. Normal forms for retarded functional differential equations 
and applications to Bogdanov-Takens singularity. J. Differential Equations 122, (1995) 
201-224. 

[10] T. Faria and L.T. Magalhaes. Normal forms for retarded functional differential equations 
with parameters and applications to Hopf bifurcation. J. Differential Equations 122, 
(1995) 181-200. 

[11] T. Faria and L.T. Magalhaes. Realization of ordinary differential equations by retarded 
functional differential equations in neighborhoods of equilibrium points. Proc. Roy. Soc. 
Edinburgh Ser. A 125, (1995) 759-776. 



13 



[12] T. Faria and L.T. Magalhaes. Restrictions on the possible flows of scalar retarded func- 
tional differential equations in neighborhoods of singularities. J. Dynam. Differential 
Equations 8, (1996) 35-70. 

[13] J.K. Hale and S.M. Verduyn Lunel. Introduction to Functional Differential Equations, 
Appl. Math. Sci., vol. 99, Springer, New York, 1993. 

[14] T. Heil, I. Fischer, W. Elsafier, B. Krauskopf, K. Green and A. Gavrielides. Delay 
dynamics of semiconductor lasers with short external cavities: Bifurcation scenarios 
and mechanisms. Phys. Rev. E 67, (2003) 066214-1-066214-11. 

[15] W. Jiang and Y. Yuan. Bogdanov-Takens singularity in Van der Pol's oscillator with 
delayed feedback. Phys. D 227, (2007) 149-161. 

[16] Y. Kuang. Delay differential equations with applications in population dynamics. Math- 
ematics in Science and Engineering, 191. Academic Press, Boston, (1993). 

[17] A. Longtin and J.G. Milton. Modeling autonomous oscillations in the human pupil light 
reflex using nonlinear delay-differential equations. Bulletin Math. Biology 51, (1989) 
605-624. 

[18] J. Sieber and B. Krauskopf. Bifurcation analysis of an inverted pendulum with delayed 
feedback control near a triple-zero eigenvalue singularity. Nonlinearity 17, (2004) 85- 
103. 

[19] E. Stone and S.A. Campbell. Stability and bifurcation analysis of a nonlinear DDE 
model for drilling. J. Nonlinear Sci. 14, (2004) 27-57. 

[20] M.J. Suarez and P.L. Schopf. A Delayed Action Oscillator for ENSO. J. Atmos. Sci. 
45, (1988), 3283-3287. 

[21] A.G. Vladimirov, D. Turaev and G. Kozyreff. Delay differential equations for mode- 
locked semiconductor lasers. Optics Letters 29, (2004) 1221-1223. 

[22] J. Wei and W. Jiang. Stability and bifurcation analysis in Van der Pol's oscillator with 
delayed feedback. J. Sound and Vibration 283, (2005) 801-819. 

[23] X. Wu and L. Wang. Zero-Hopf bifurcation for van der Pol's oscillator with delayed 
feedback. J. Comput. Appl. Math. 235, (2011) 2586-2602. 

[24] P. Yu and Y. Yuan. The simplest normal forms associated with a triple zero eigenvalue 
of indices one and two. Nonlin. Anal. 47, (2001) 1105-1116. 



14 



[25] C. Zhang and J. Wei. Stability and bifurcation analysis in a kind of business cycle model 
with delay. Chaos Solitons Fractals 22, (2004) 883-896. 



15 



